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Abstract 

We study the dynamics of a suspension of magnetic nanoparticles. Their relaxation times are 
strongly size-dependent. The dominant mode of relaxation is also governed by the size of the par- 
ticles. As a result the dynamics is greatly altered due to polydispersity in the sample. We study 
the effect of polydispersity on the response functions. These exhibit significant changes as the 
parameters characterizing polydispersity are varied. We also provide a procedure to extract the 
particle size distribution in a polydisperse sample using Cole-Cole plots. Further the presence of 
attractive interactions causes aggregation of particles leading to the formation of clusters. Repul- 
sive interactions along with thermal disorder not only hinder aggregation, but also introduce the 
possibility of removal of particles or "fragmentation" from clusters. The competing mechanisms of 
aggregation and fragmentation yield a distribution of cluster sizes in the steady-state. We attempt 
to understand the formation of clusters and their distributions using a model incorporating the 
phenomena of aggregation and fragmentation. Scaling forms for quantities of interest have been 
obtained. Finally we compare our numerical results with experimental data. These comparisons 
are satisfactory. 

PACS numbers: 47.57.E-,47.65.Cb,47.57.eb,89.75Da 
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I. INTRODUCTION 



Single domain magnetic nanoparticles (MNP) and their colloidal suspensions have at- 
tracted a lot of attention in the recent years H-Qj- The growing interest is due to a variety 
of technological applications associated with them. These range from mechanical and ther- 
mal applications involving their usage as sealants, lubricants and coolants to challenging 
applications in medicine for the purpose of magnetic resonance imaging, targeted drug de- 
livery and biomarkers and biosensors to name a few. The main reason behind their wide 
applicability is the ease with which they can be detected and manipulated by the applica- 
tion of an external magnetic field. Their response times are strongly size-dependent, thus 
introducing the possibility of synthesizing particles to yield application tailored response 
times. 

Most practical applications require appropriate surfactant coatings to prevent agglomera- 
tion and sedimentation of magnetic nanoparticles. Many biological and medical applications 
require nanoparticles with biologically relevant coatings in order to use them as probes and 
carriers j^, 5|. Such coatings enhance the particle size, although the magnetic volume re- 
mains unchanged. The Neel and Brownian relaxation times which characterize the dynamics 
of the suspension depend not only upon the constituting material, but also the magnetic 
volume and the enhanced volume due to surfactant coating. The interplay of theseparame- 
ters governs whether one or both relaxation times contribute to the dynamics In 
the present paper, we systematically analyze the effect of the above parameters on the relax- 
ation times. We also study the ac susceptibility the most commonly studied response 
function, in various regimes characterized by the relative dominance of the two relaxation 
times. This understanding makes it possible to estimate relaxation times as well as related 
parameters from the measurement of xi u ) i n t ne laboratory. 

All experimental samples have a distribution of particle sizes and are referred to as 
polydisperse. The particle size distributions in experimental samples, usually obtained by a 
tunnelling electron microscopy (TEM) analysis, are found to have a log- normal form l|, |8[. 
Hence we incorporate the effect of polydispersity in our calculations of As the Neel 

and Brownian relaxation times have a strong dependence on particle size, polydispersity 
leads to a considerable broadening and sometimes an additional peak in We have 

also worked out a procedure to obtain particle size distributions from the susceptibility data 
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via Cole-Cole plots when they are unavailable 10, 111. 

Q 

The above approach assumes the single particle model [7J applicable to dilute suspen- 
sions in which magnetic interactions amongst particles may be ignored. Many biological 
applications are particularly benefitted by this approach. For instance diagnostic tools such 
as magnetic resonance imaging rely on transport and manipulation of individual magnetic 
particles bound by fluorescent dyes through blood vessels [l2|. However many mechanical 
and thermal applications rely on viscosity effects arising due to a higher density of MNP in 
the suspensions. When present in sufficient concentration, clustering and chaining of MNP 
is rather common as observed by electron microscopy or dynamic light scattering studies 
[13I ]- 15| . This behavior is undesirable in many applications using magnetic fluids as sealants, 
coolants, lubricants, printing inks, etc. where invariance of the magnetic and fluid properties 
are paramount [1]. On the other hand magnetic domain detection , op tical shutters, tagging 
of surfaces and other entities benefit from clustering of MNP [l|, [lg]. It is hence useful to 
understand the mechanisms responsible for clustering and the dependence of cluster-size dis- 
tributions and average cluster sizes on the parameters describing the experimental sample. 
With this knowledge, it may be possible to synthesize application tailored suspensions. 

Aggregation results due to a variety of interaction energies which come into play in mag- 
netic suspension jl[ I2I, Q . Typically these are dipolar, van der Waal and steric interactions. 
Of these the former two are attractive while the later is repulsive. Their relative strengths 
are governed by the magnetic compound, the magnetic volume, the thickness of the surfac- 
tant coating as well as the properties of the suspending liquid. Additionally, temperature 
acts as a disordering agent which can not only hinder aggregation but also introduce the 
possibility of detachment of the basic unit(s) from the parent cluster. Thus we expect the 
cluster-size distribution and average cluster size to be governed by the interplay of the at- 
tractive and repulsive forces amongst the interacting particles and the ambient temperature 
of the suspension. 

Earlier models addressing the issue of clustering in MNP suspensions have treated ag- 
gregation as an irreversible mechanism. These models obtain a power law growth for the 
average cluster size in the initial stages of aggregation while the steady-state is characterized 
by a single cluster comprising of all the particles in the suspension [^-[^j]. The later is 
often referred to as an infinite cluster. We further add fragmentation in the model to counter 
irreversible aggregation, due to the presence of repulsive interactions and temperature. The 
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two competing mechanisms make cluster formation a self-limiting process resulting in a 
distribution of clusters of varying sizes with the average cluster size governed by the ra- 
tio R of the aggregation and fragmentation rates. We obtain the steady-state cluster-size 
distributions and the evolution of the mean cluster size to its steady-state value for differ- 
ent values of R. Both these functions exhibit scaling. We have compared our results with 
sets of experimental data to lend credence to our interpretations and observations of the 
aggregation-fragmentation model. These comparisons are much to our satisfaction. 

The paper is organized as follows. Section [III deals with relaxation mechanisms and re- 
sponse functions of dilute MNP suspensions using the single particle model. In Section III A[ 
we introduce the Neel and Brownian relaxation times and identify regimes where either or 
both relaxation times contribute. Polydispersity and its characterization is introduced in 
Section [TTBI The calculation of the ac susceptibility f° r monodisperse and polydisperse 
samples in the different regimes described above is presented in Section Hi CI In Section Hi D\ 
we provide a procedure to obtain particle-size distributions (for polydisperse samples) using 
Cole-Cole plots obtained from Section [TTT1 deals with aspects of clustering occuring 

in samples which are no longer governed by the single particle approximation. Interactions 
prevalent in samples exhibiting these phenomena are introduced in Section IIII Al We in- 
troduce the aggregation-fragmentation model to describe clustering in Section IIII BL The 
numerical results are presented in Section nil CI Their comparisons with corresponding mea- 
surements in a variety of experimental systems are presented in Section IIII Dl Finally we 
conclude this paper with a summary and discussion in Section [TV] 

II. RELAXATION MECHANISMS AND RESPONSE FUNCTIONS OF DILUTE 
MNP SUSPENSIONS 

We now study the relaxation properties and the response functions in dilute MNP suspen- 
sions where inter-particle interactions can be ignored. In this regime the particles essentially 
behave as independent, single- domain, super paramagnetic entities. 
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A. Neel and Brownian relaxation times 



A ferromagnetic sample such as iron for instance, has no net magnetization even below 
the Curie temperature. This is so because the sample comprises of domains, each having 
spontaneous magnetization pointing in a different direction. However if the size of the body 
is reduced, there comes a point beyond which the magnetostatic energy and the energy for 
orming domain walls compete in such a way that a single domain state becomes preferable 



231 ] . The direction of magnetization of the single domain particle does not remain fixed in 
time though, but undergoes fluctuations or "relaxations" as the magnetic moment rotates 
between the crystallographic anisotropy axes. As a result the time averaged magnetization 
is still zero and the particle is still paramagnetic. It is called "super paramagnetic" because 
each particle has a giant magnetic moment arising due to a large number (~ 10 5 ) of individual 
atomic moments. We now describe the mechanisms of relaxation of these single domain super 
paramagnetic particles suspended in a liquid. 



The magnetization vector of a single domain particle is given by 23]: 



M = VM h, (1) 

where V = 4nr^ /3 is the magnetic volume of a particle with radius r c usually reffered to as 
the core radius, M Q is the saturation magnetization and n a unit vector in the direction of 
the magnetization. In the case of uniaxial anisotropy (in the z-direction say), the magnetic 
energy is given by 

E = VK sin 2 6, (2) 

where K is the effective magnetic anisotropy constant and 9 is the angle between the z- 
axis and h. Minimum energy occurs at 9 = and it defining two equilibrium orientations 
corresponding to magnetizations +VM D and —VMo. If thermal fluctuations are strong 
enough, magnetic moment reversal takes place within the particle by overcoming the energy 
barrier (of height VK). This reversal or switching time is called Neel relaxation time and 
is given by 

r N = r e VK ^ T , (3) 

where r G is related to the inverse of the attempt frequency of magnetic reversal. 

There is another mechanism by which the magnetic moment of a super paramagnetic 
particle suspended in a fluid can relax. This mechanism of relaxation can be due to the 



physical rotation of the particle within the fluid. It is referred to as Brownian rotational 
motion as it occurs due to the thermal fluctuations in the suspended medium. The Brownian 
relaxation time is given by {^J 

= -J^r> (4) 

where r\ is the dynamic viscosity of carrier liquid and is the hydrodynamic radius defined 
as the sum of the core radius r c of the MNP and the surfactant coating 5 over it. 

As can be seen from Eqs.(j3]) and (111) , both Neel and Brownian relaxation times are highly 
sensitive to the particle size. While increases exponentially, tb grows linearly with the 
particle dimension. It is customary to define an effective relaxation time as follows {7, 8]: 

1 1 1 

— = 1 or 

T e T N T B 

TN + T~B 

T e = ~ ~ • ( 5 ) 
TnTb 

Thus it is possible to tailor time scales by an appropriate choice of parameters, particularly 
K, r c and r^. With this in mind, we have systematically studied the effect of r c and on 
ttv and tb and consequently r e for maghemite {Fe^Oi) particles used most commonly in 
making magnetic fluids. 

In Table I, we summarize our evaluations of Tn, Tb and r e for particles with varying 
magnetic core radius r c and the thickness of the surfactant coating 5. In most experiments, 
the later is usually in the range of 2-6 nm. As is observed in Table I, tb is practically 
unaffected by 5. For small particles, rjy <C tb which results in r e « rjy. The relaxation then 
takes place by rotation of the magnetic moment inside the particle. For large particles on 
the other hand, tb <C t~n- Consequently r e w t b and the relaxation is due to a physical 
rotation of the particle in the suspension. Thus the choice of the relaxation mode is primarily 
governed by the particle size. For critical particle sizes, often called the cross-over radius r* 
(~ 8 nm in Table I), it is found that both mechanisms contribute to the relaxation of the 
suspended particle. 

We have studied the effect of temperature T, the surfactant coating 5 and the anisotropy 
constant K on the cross-over radius r*. Most applications require an operating temperature 
in the range of 270 to 320 K. We find that r* does not change perceptibly in this range. 
Further, as seen from Table I the variation of 5 does not significatly alter r*. The anisotropy 
constant K on the other hand, leads to a substantial change in the corresponding value of 
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the cross-over radius. In Figure 1, we plot r* as a function of the anisotropy constant K 
for a few frequently used biological and commercial spherical, magnetic nanoparticles. In 
all these evaluations, the temperature is assumed to be 300 K and the surfactant coating 5 
has been taken to be 2 nm. The figure indicates that the value of the cross-over radius r* is 
smaller for larger values of anisotropy K. Thus amongst the three parameters of relevance, 
the anisotropy constant K affects the cross-over radius r* the most. 



B. Polydispersity 

Monodisperse samples are an idealization. A variation in the particle size is inherent 



in all experimental samples. TEM studies o: 



distribution for the variation of particle sizes 



several samples have revealed a log-normal 



261 ] . Thus the probability density P(r c )d 



of having particles within radius r c and r c + dr c can be written as: 

P{r c )dr c = —J- exp [- In 2 (r c /f c )/ (2 In 2 a)] dr c , (6) 

V 2tt in a 

where f c and a are the mean and variance of the distribution. Due to the strong dependence 
of both Neel and Brownian relaxation times on particle size, it is evident that a distribution 
of relaxation times will be obtained if the suspended particles have a distribution of sizes. 
The presence of polydispersity leads to a significant change in the behavior of the response 
function of the sample as we shall see in the following subsection. 



C. AC susceptibility measurements 

The complex susceptibility °f a suspension of monodisperse MNP in the linear 

response regime has the Debye form given by [231 ]: 

XM = t: X ° v (7) 

(1 - LUT e ) 

where Xo = x( w = 0) = NV 2 M^/ksT is the static susceptibility of the sample comprising 
of N monodisperse particles of volume V with a saturation magnetization M Q . The effective 
relaxation time r e is defined by Eq. flS}. An experimental time scale is provided by the 
inverse of the frequency of the applied oscillatory field. If ur e 1, the response is given by 
the static susceptibility Xo bearing the characteristic T _1 dependence. On the other hand, 
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if the two time scales are comparable, i.e. when wr e m 1, the response is marked by strong 
frequency dependent effects. Separating the real and the imaginary parts of susceptibility 
yields: 

x'M = Xo : , 1 , , and (8) 

1 + U) 2 T* 

x» = xo -^r- r (9) 

The imaginary part of the susceptibility, x"(cu) governs dissipation in the system. When 
plotted as a function of the frequency, it exhibits a symmetric peak around u = t~ x . Al- 
ternatively, the peak frequency can provide information about the radius r c as well as the 
hydrodynamic radius of the particle as will be discussed shortly. 

The susceptibility response gets substantially altered in the presence of polydispersity. 
For polydisperse samples, needs to be averaged over the particle size distribution 

P(r c ). Thus 

XV) = Xo I dr c P(r c )— — \ and (10) 

J l + u 2 r 2 (r c ) 

X » = Xo / d r c P i r < .) T ^>_. (11) 

In Figure 2, we plot Eq. (Q (open circles) for monodisperse samples corresponding to three 
different values of r c = 4, 8 and 12 nm. The chosen values correspond to (a) r c < r*, (b) 
r c m r* and (c) r c > r* respectively for maghemite particles used to generate data of Table I. 
All cases exhibit the characteristic Debye form. To understand the effect of polydispersity, 
we also plot Eq. (II ip (filled circles) in the same figure. The distribution P(r c ) vs. r c used for 
;he evaluations were obtained from a TEM analysis of maghemite samples used in reference 
8J. These distributions had a log-normal form with a variance a ~ 0.35. The mean value 
r c is the particle size of the corresponding monodisperse evaluation. 

As seen in Figure 2, a broadening of the response function is observed in all cases after 
the inclusion of polydispersity. The response no longer has the symmetric Debye form. In 
Figure 2a, the peak frequency provides information about the mean core radius f c . The 
width of the log-normal distribution often gives rise to a small Brownian peak although f c < 
r* and the particles predominantly exhibit Neel relaxation A pronounced two-peak response 
is obtained in Figure 2b since r c ~ r*. Both Neel and Brownian relaxation contribute 
in this regime. The frequencies corresponding to the peak values yield information about 
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the average values of core and the hydrodynamic radii f c and f~h of the particles. Finally 
when r c > r*, the Brownian relaxation dominates and the frequency corresponding to the 
peak provides information regarding the hydrodynamic radius of the particle. Comparing 
Figures 2a and 2c, it is clear that polydispersity affects Neel relaxation more significantly 
than Brownian relaxation as expected (cf. Eqs. ([3D and 



HQ: 



*H = ^T^-T^, (12) 



D. Particle size distributions from Cole-Cole plots 

The information on distribution of relaxation times can be obtained from susceptibility 
measurements by empirical models. The most frequently used model for obtaining the 
distribution of relaxation times is the Cole-Cole model. Based on this model, the dynamic 
susceptibility x{ u ) with multiple relaxation times is given by 

Xo 

1 + (iOJT c )' 

where r c is the central relaxation time about which all the other relaxation times are dis- 
tributed and a is a fitting parameter with limits 0< a <1. The equation reduces to the 
Debye equation for a=0. As the deviation from the single relaxation time model becomes 
greater, a — >1. Separating the real and imaginary parts of Eq. (fl2|) . we obtain: 

X'H = ^ U - w SinMl ~^ and (13) 

Ay ' 2 V cosh(l - a)s + cos(avr/2) ) v ; 

y"(u) = ^ cos(a7r/2) (U) 

XK ' 2 cosh(l-a)s + sin(a7r/2)' 1 ; 

where s = logo;r c . 

Cole and Cole proposed a method of graphically representing the effects of multiple 
relaxation times. The method consists of plotting x"( w ) f° r a certain frequency against 
x'(w) at the same frequency. These are called Cole-Cole plots. When a = 0, the Cole-Cole 
plot is a semi-circle. When a > 0, the Cole-Cole plot is still a semi-circular arc similar 
to a widened Debye curve, but in this case the centre lies below the horizontal axis. The 
plot is symmetrical about the vertical line passing through the point x'{ u ) = Xo/2 when 
x"( w ) is maximum at a frequency u> = t~ x . The parameter a can be determined from the 
Cole-Cole plot by a graphical construction. We do not reiterate this rather well established 
procedure here, but refer the reader to reference [10] for it. It may also be determined by 
best fits of Eq. ffl3|) and Eq. f JT4|) to the experimental data. It should be noted that the 



symmetrical distribution of relaxation times is a consequence of a particle size distribution 
which is symmetric about the mean value f c . The Cole-Cole model assumes a gaussian form 
for P(r c ) vs. r c . The parameter a is hence related to the spread of the gaussian distribution 
with a = yielding a delta function or a monodisperse particle size distribution. 

The response functions in typical MNP suspensions, as observed from Figure 2, are 
asymmetrical primarily due to the log-normal particle size distributions. The Cole- Cole 
model is thus not the most suitable for their description. Experimental systems exhibiting 
asymmetric response functions can be conveniently represented by an expression due to Cole 



and Davidson. The ac susceptibility in this model is given by 



(1 + IUJT C )P 

where < /3 < 1 and is related to the deviation from the single relaxation time model. The 
Debye form of Eq. ([7]) is recovered when /3 = 1. Separating the real and imaginary parts of 
Eq. (fl5l) results in: 

x'(^) = Xo cos <j/ 'cos /3(f) and (16) 
= Xo cos <$P sin p<f>, (17) 

where = arctano;r c . Eqs. (TIB"]) and ( TTT1) give rise to a "skewed arc" in the Cole-Cole plots. 
The parameter can be obtained by a graphical construction or by fitting the experimental 
data corresponding to %!(co) an d x"^) with Eqs. ffT6|) -f lT7|) . 

In Figure 3 we show the Cole-Cole plots for the same polydisperse samples that were used 
to obtain data of Figure 2. The cases (a), (b) and (c) correspond to r c < r*, r c r* and 
f c > r* respectively. In all the three cases, the centre of the arc lies below the horizontal axis 
due to multiple relaxation times resulting from polydispersity. Further, the asymmetry in 
the particle size distributions is reflected in the asymmetry of the arcs. The Cole-Cole plots 
are distinct in each of the regimes. In particular, the cross-over regime is characterized by a 
two- humped form signifying a comparable contribution of both Neel and Brownian relaxation 
to the response function. Thus significant information regarding the experimental sample 
can be inferred from the Cole-Cole plots. 

We now provide a simple procedure for evaluating the particle-size distribution P(r c ) 
vs. r c in the experimental sample when the response function exhibits multiple relaxation 
times. Assuming a log- normal form, P(r c ) vs. r c is characterized by its mean r c and 
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variance a signifying the spread in the particle-size distributions. While the evaluation of f c 
is straightforward, the variance a needs to be estimated. Recalling that multiple relaxation 
times are a consequence of polydispersity, it is imperative to connect j3 with a. In order to 
find this relationship, we have gone through the following sequence of steps. Firstly x'i 00 ) 
and were evaluated using Eqs. (fT0|) and (TTTT) for a chosen value of r c and a. The 

corresponding values of r c and (3 were then obtained by fitting the above susceptibility data 
with Eqs. ffTB"]) and f|T7|) of the Cole-Davidson model. This procedure was then repeated for 
a number of o values in the range 0.1 to 0.5, keeping r c constant. Larger values of a were 
not considered as they resulted in extremely broad log-normal distributions which are not 
of relevance in experiments. The above evaluation was then performed for three different 
values of K. The results obtained are plotted in Figure 4 on a semi- logarithmic scale. The 
solid lines follow the equation: 

a = A(K)e- B P. (18) 

The parameters A (K) and B in each case are obtained from the best fits of the above 
equation to the evaluated data. We find that B = 2.45±0.05 and is independent of K. We 
also find that Eq. (1181) and the parameter B is unaffected by r c , which was varied from 4 to 
12 nm. 

We test the procedure for extracting particle-size distribution from the susceptibility data 
on maghemite nanoparticle suspensions provided in reference 8] . From the Cole-Cole plots 
of this data, we evaluated the best-fit values of r c and 0. These were found to be 2 x 10 -8 
s and 0.44 respectively. The average core radius r c , evaluated from r c using Eq. ([3]), was 
found to be 4.2 nm. The value of a characterizing the spread of the log-normal distribution 
was then obtained using Eq. ffTB"]) and was found to be 0.36. The reconstructed particle- 
size distribution is shown in Figure 5. The particle-size distribution obtained using TEM 
measurements has also been provided in the figure for comparisons. In general, we find that 
this procedure reproduces the original distributions rather well when the response function 
has one dominant peak, be it in the Neel or the Brownian relaxation regime. 

III. CLUSTERING IN MNP SUSPENSIONS 

The properties of the magnetic fluid are greatly affected by the aggregation of particles 
(inspite of surfactant coating) due to the presence of attractive and repulsive interactions. 
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It is hence instructive to consider the different interaction energies in MNP suspensions and 
understand their interplay which leads to the formation of aggregates or clusters of different 
sizes. 



A. Interactions in MNP suspensions 

The primary interaction energies in these systems are enumerated below 

1. Dipolar interaction: 

As each particle is magnetized, the attractive dipole-dipole interaction could in princi- 
ple force particle agglomeration. The dipolar interaction energy between two magnetic 
nanoparticles, each having a magnetic moment fi is given by: 

4 7T \ S° S 6 J 

where s is the center-to-center separation between the two nanoparticles and the per- 
meability of free space fi Q = Air x 10~ 7 Hm~ x . As the magnetic properties of the 
particles are affected by temperature, the dipolar interaction defined above is also 
temperature dependent. 

2. van der Waal's interaction: 

It arises spontaneously between neutral particles because of the fluctuating electric 
dipolar forces and is attractive in nature. Hamaker calculated this interaction for 
identical spheres separated by a surface-to-surface distance / to be: 

In the above equation, the Hamaker constant A = 10~ 19 Nm. It should be noted that 
the / _1 dependence in the above equation indicates that infinite energy is required to 
separate a particle pair while a finite energy is required in its formation. Therefore 
agglomeration of particles will occur as long as the Hamaker constant has a finite 
value. As a result van der Waal's interaction unlike dipolar interaction is unaffected 
by temperature. 

3. Steric interaction: 

The steric energy comes into play due to the presence of long chained surfactant 
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molecules coating the particles. This mechanism prevents the particles from approach- 
ing very close to one another, thus preventing the van der Waal's attraction to come 
into play. This repulsive interaction originates due to a compression of the surface- 
adsorbed surfactant molecules when the inter particle separation is smaller than two 
times the thickness of the surfactant layer. For a sufficient density of the surfactant 
layer, the repulsion can grow large enough to avoid the contact between the magnetic 
particles. This repulsive energy for spherical particles has been calculated by Papell 
and is given by the following form: 

E. = ^-(2S-i) 1 + (21) 
1.325 V 2/ V 2 / 

In the above form, the Young's modulus E was assumed to be 10 6 and d c = 2r c is the 
diameter of the magnetic core of the particle. 



Apart from the above interaction energies, the thermal energy responsible for the Brow- 
nian motion of the suspended particles also plays an important role in the aggregation 
dynamics. An appropriate value of temperature (which depends on the size of the inter- 
acting particles) is capable of preventing aggregation. In principle, this condition may be 
expressed as a comparison between the attractive energy with the thermal energy of the two 
particles which yields a numerical value of the core radius r c at which this could happen. 
We expect the thermal degrees of freedom, in conjunction with steric repulsion, to not only 
hamper the aggregation process but also remove particles from the parent cluster leading to 
its fragmentation. 

At this juncture we define a parameter Q as the ratio of energies leading to aggregation 
to those which lead to fragmentation follows: 

n = Ed + K . (22) 

E s + k B T 1 ; 

It is expected that the distribution of cluster sizes and the mean cluster size will be governed 
by the ratio Q. In the following subsection, we present a model which incorporates the 
competing mechanisms of aggregation and fragmentation. 
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B. The aggregation-fragmentation model 



The formulation introduced by Smoluchowski is especially useful to model a suspension 



of magnetic nanoparticles 



19j. To begin with, we assume that the suspension contains N 



identical, single particles executing Brownian motion. The later leads to aggregation of two 
particles if they come within an appropriate range of one another and the net interaction 
between them is attractive. The resulting cluster of size two also executes Brownian mo- 
tion, but with a reduced diffusion rate till it encounters a particle or a cluster of particles. 
The process goes on and eventually a single large cluster comprising of all the N particles 
is formed. Such a cluster is usually reffered to as an infinite aggregate. In most useful 
suspensions however there is a distribution of clusters of varying sizes. As discussed in the 
preceding subsection, the combined effects of thermal energy and inter-particle repulsion can 
introduce fragmentation in the cluster dynamics thereby preventing the formation of large 
aggregates. We thus include this additional mechanism in the rate equations which describe 
the evolution of clusters. 

Let c(k,t) denote the number of clusters containing k particles at time t. The time 
evolution of c(k,t) is governed by the following rate equations: 



= ~ K ijC {i,t)c{j,t) - c(k,t)jrK k]C (j,t) 

i+j=k j=l 

oo 

+fwc(k + l,t)-f k c(k,t) + 6 ktl ^T l f j c(j,t), k>l. (23) 

i=i 

In the above equation the Kij and ft are the aggregation and fragmentation kernels respec- 
tively whose forms are specified below. The aggregation kernel describes the coalescence of 
a cluster containing i particles with another containing j particles to yield a larger aggregate 
comprising of k = i + j particles. It is assumed to have a mass- dependent form defined by 
= D{i~^ + to take into account the reduced mobility of large clusters. A value of 
[i = implies a mass- independent mobility, i.e., clusters of all sizes diffuse with the same 
ease. A non-zero values of /i results in a slower mobility of larger clusters and consequently 
a slower growth rate of clusters. In the limit of fi = oo, only monomers are mobile. The 
choice of [i is dictated by the experimental parameters as we shall see in a short while. The 
fragmentation kernel fk describes the loss of a particle from the parent cluster and can also 
be assumed to have a mass-dependent form defined by fk— oj k v . In our studies however, 
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we have set v = to reduce the number of parameters in the model. The parameters D and 
u define the relative strength of the aggregation and fragmentation processes. 

The first and the third term in Eq. ()23p are referred to as the gain terms which result 
in the formation of clusters of size k. The former describes aggregation of two clusters to 
yield a cluster containing k particles while the later describes the generation of a cluster 
comprising of k particles due to fragmentation of a particle from a cluster of size k + 1. The 
second and the fourth term on the other hand are referred to as the loss terms which describe 
processes leading to loss of clusters of size k. This could be due to aggregation of a cluster 
of size k with another or the its fragmentation. The fifth term in the equation describes the 
generation of single particles from clusters due to the process of fragmentation. It is easy to 
check that the following sum rule is satisfied: 



0_ 

Of 



J2 kc (k,t) = 0, or J2kc(k 1 t) = N, (24) 

,fc=l / k=l 

as required by conservation conditions. 

The condition dc(k, t)/dt = describes the steady-state which is of interest to us. In the 
absence of the fragmentation term, Eq. (123]) reduces to the Smoluchowski equation describing 



coagulation phenomena 19|. There have been a few studies of this model, both analytical 
and numerical, to predict scaling forms associated with cluster growth and cluster-size dis- 
tributions {20)]- 22], 28]- 31 1 - The growth of clusters with time is a power governed by the 
relation < k(t) > ~ t z . Choosing ji — 1, the value of z was estimated to be ~ 0.5. The 
steady state in this model was found to be an infinite (single) aggregate comprising of all 
the N particles. Note that in the absence of fragmentation, the constant D can be absorbed 
by redefining t as Dt in Eq. ( I23p . making the later independent of the aggregation rate. 

To mimic those physical situations which do not have an infinite-aggregate as a steady- 
state but rather have a distribution of clusters of varying sizes, it is essential to include 
fragmentation as a competiting mechanism to aggregation. Of relevance in the context of 
the present study is model incorporating mass- independent aggregation and mass-dependent 



evaporation processes reported in reference 29[. The competition between aggregation and 



evaporation leads to several asymptotic outcomes of the steady-state solution of this model. 
For instance if evaporation dominates over aggregation, the steady-state cluster size distri- 
bution c(k) vs. k decays exponentially. On the other hand for a critical evaporation rate, 
the distribution decays as k~ 5 ^ 2 . 
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It is convenient to rewrite Eq. (!23|) in terms of P(k, t), the probability of having cluster of 
size k at time t. Defining the later as P(k, t) = c(k, t)/ J2T=i c (^' 0? ^s evolution is governed 
by the following equation: 

= K a p ^ *) p o> t) - p(k, t) jr K kj p(j, t) 

i+j=k j=l 

+wP(k + l,t)-wP(k,t), k>l, (25) 

^M = -p(u)f;^?( jl t)+ w f;p( j ,t) 1 k = i. (26) 



The above equations also satisfy the sum rule 

' oc \ oo 

J2 P (k,t) = 0, or £p(*,t) = l, (27) 



9t 



v fc=i / fc=i 



as required. 



C. Numerical results 



We now solve the the set of equations defined by Eqs. (!25l) - (l26l) numerically to obtain 
the steady-state cluster size distribution P(k) vs. k. It is useful to define the ratio R 
characterizing the relative strengths the aggregation and fragmentation mechanisms: 

R = —. 28 

Identifying the physical origin of D in the attractive interactions between clusters and that of 
u in the disordering agents (such as repulsive interactions and temperature), we can expect 
R to have the same qualitative effect as the ratio Q defined in Eq. (1221) on steady-state 
cluster size distributions. We use this correspondence to bring contact between numerical 
results and experimental observations on cluster formation in MNP suspensions. 

In Figure 6a, we look at the variation of the mean cluster size < k(R,t) > as a function 
of t for different values of R on a double logarithmic scale. The parameter \x was chosen to 
be 2.0. After an initial growth period obeying a power law, the cluster size attains a steady- 
state value of < k(R) > due to the competing mechanisms of aggregation and fragmentation. 
As expected, < k(R) > increases with increasing values of R. In fact we find that < k(R) > 
~ R a with a = 0.85 ±0.02 as can be seen in the inset. Further, in the diffusion domiated 
regime (R > 0.5) the data in Figure 6a can be scaled by replotting < k(R,t) > / < k(R) > 
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vs. t/t s where t s is the time taken to attain the steady-state. The initial cluster growth is 
of the form < k(t) >~ t z . We find that the growth exponent z = 0.38 ± 0.02 . The scaled 
data is shown in Figure 6b on a double logarithmic scale. The best fit line with a slope of 
0.38 is also indicated. The value of the growth exponent also depends upon the value of 
/i. For n = 1.0, our simulations yield z = 0.78 ± 0.02. The corresponding scaled data and 
the power law fit is also shown in the inset of Figure 6b. The faster growth of clusters is a 
consequence of increased mobility due to a lower value of /i. 

Next in Figure 7a we plot the steady-state distribution P{k) vs. k corresponding to fi = 
2 for the same set of R values considered above. We find that the tails of the distributions 
fit well to a power law in k. In Figure 7b, we plot < k(R) > P{k) vs. kj < k{R) > where 
< k(R) > is the steady-state average cluster size for the corresponding value of R. These 
plots indicate that the distribution functions P{k) corresponding to different values of R 
obey scaling in the diffusion-dominated regime. The scaling relation can be summarized in 
the following equation: 

P « = TltW fei^) ■ <*<*»-*■■ P 9 ) 
D. Comparisons with experimental data 

We now compare our simulation results with experimental data on average cluster sizes 
in a variety of magnetic nanoparticle suspensions. Eberbeck, et al studied the aggregation 
of various magnetic nanoparticles in a variety of media such as water, phosphate buffered 
saline, calf serum, bovine serum and human serum 17], |32|. The experimental measurements 



were performed at room temperature. In all cases, formation of dimers and trimers were 
reported by the authors. To make comparisons with the simulation results of our model, 
we first calculate the ratio Q defined in Eq. (|22p for each of the samples. These along with 
the sample specifications are tabulated in Table II. As can be seen in column 4, Q is in the 
range 0.35 to 0.70. Referring to the plot of < k(R) > vs. R for fi = 2 in the inset of Figure 
6a, a value of R in the above range yields dimers and trimers in the steady-state. 

Next we look at experiments which report the formation of long chains of magnetic 
nanoparticles in the presence of an applied magnetic field H. The particles acquire dipole 
moment \x given by: 

4 

P> = ^rl/ioxH, (30) 
17 



where r c is the core radius of the magnetic particle and x is the magnetic susceptibility. The 
interaction energy between two particles with aligned, identical dipole moments is given by: 

tti „s 1 — 3cos 2 6 , . 

u m = ■££—?—• (31) 

where s is the centre to centre separation between the two particles and 9 is the angle 
between the applied field and the line joining the centres of the spherical particles. The 
interaction is thus attractive when the dipoles are head-to-tail and repulsive when they are 
side-by-side. For these samples, the dominant interaction is dipolar. Thus the parameter Q 
defined in Eq. (122]) reduces to 

(32) 

where we have assumed the dipolar interaction to be at its maximum, when the particles are 
aligned with the external field. In the presence of an external field, if the resulting magnetic 
dipole-dipole interaction between particles exceeds the thermal energy, then the particles 
arrange themselves into chains which are column-like for low volume fractions and worm-like 



structures for high-volume fraction 



13 



14|. Many theoretical studies on the above 



systems 



have been based on Smoluchowski kinetic equations for irreversible aggregation 20l. 121]. 

A set of experimental measurements which we find especially relevant in the context of the 
aggregation- fragmentation model are reported in reference [13]. In this paper, the authors 
have studied aggregation dynamics in very dilute emulsions of ferrofluid droplets in water. 
The ferrofluid droplets were small Fe^O^ grains in kerosene coated with a surfactant to 
prevent agglomeration. The data on average chain lengths as a function of time for different 
volume fractions and applied fields was obtained using dynamic light scattering experiments. 
We reproduce a scaled form of this data in Figure 8. On the x-axis, we plot t/t s where t s 
is the time taken to reach the steady-state value of the average cluster. The y-axis has 
been scaled by the steady-state value of the average chain length. The unsealed data is also 
provided in the inset for reference. As can be observed in Figure 8, the scaled form (as well 
as the unsealed form) of the experimental data bears a qualitative resemblance to Figure 6b 
which results from the numerical solution of the aggregation-fragmentation model defined 
by Eqs. ( 1251) and ( )26|) . The slow initial growth, the power law form at intermediate times 
and saturation to a steady-state value as time progresses are bourne by both sets of data. 

On the quantitative side, we find that the experimental data is well represented by a 
growth exponent z = 0.78±0.02 obtained with fi = 1 in our simulations. The corresponding 
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line is depicted in figure 8. It is pertinent to recall here that the dipolar interaction becomes 
dominant in the presence of the applied magnetic field thereby increasing the diffusivity 
of the clusters. In the aggregation-fragmentation model, the later can be achieved by an 
increase in the value D as well as a decrease in the value of fi. Consequently, we find that 
our data on < k(R,t) > vs. t for a value of \i = 1 rather than \x = 2 results in accurate 
comparisons with experimental data. The initial growth does not conform to the predicted 
value of z — 0.5 (also obtained with \i — 1) of the irreversible aggregation model. 

A similar set of experimental results were obtained in reference 14| where the power law 
prediction for cluster growth was checked by performing a large number of experiments on 
aqueous solutions of super paramagnetic polystyrene beads having a uniform distribution 
of Fe^O^ particles. The data was obtained for five values of the volume fraction in the 
presence of low field strengths and higher field strength using optical microscopy. The 
exponent z describing the cluster growth in these studies was found to be much larger than 
the predicted value of 0.5 especially for small volume fractions and low field strengths. These 
observations further reiterate the appropriateness of introducing fragmentation along with 
aggregation especially in the above limits where the thermal energy plays a significant role 
in the dynamics. 



IV. SUMMARY AND DISCUSSION 



We conclude this paper with a summary of results and discussion presented here. Our 
main interest was to understand the factors governing the dynamical response of suspensions 
of single-domain magnetic nanoparticles. Such an understanding introduces the possibility 
of synthesizing particles with application tailored response times. The effect of sample pa- 
rameters on the Neel and Brownian relaxation times which characterize the response was 
studied. Amongst all the parameters of relevance, the anisotropy of the constituting mate- 
rial and the particle size alter the relaxation time most significantly. Infact the dominant 
relaxation time is also decided primarily by the particle size. We also studied how these 
parameters affect the ac susceptibility which is the most commonly studied response 
function in the laboratory. This understanding proves to be useful in estimating relaxation 
times as well as sample parameters from the measurement of x{ u ) i n the laboratory. 

We have then studied the effect of polydispersity, an inherent feature of all samples on the 
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response characteristics. These exhibit significant changes due to the strong dependence of 
relaxation times on particle sizes. The primary effect is the broadening of response functions 
and in some cases the later exhibits a two peaked structure. We have also worked out a 
procedure to obtain the particle-size distribution from x{ u ) using Cole-Cole plots and the 
analysis of Cole and Davidson. This provides an alternative approach to TEM analysis 
which is usually employed to obtain particle-size distributions. 

The above studies assumed a single-particle model applicable to dilute suspensions. How- 
ever in many cases the inter-particle interactions cannot be ignored. Clustering of particles is 
inevitable in such suspensions. While the formation of clusters is undesirable in some appli- 
cations, it is beneficial in many others. Hence we have tried to understand the mechanisms 
responsible for clustering and the experimental parameters which govern the properties of 
cluster-size distribution and the average cluster size. This knowledge is useful for the syn- 
thesis of application tailored suspensions. 

A model incorporating the phenomena of aggregation and fragmentation was used to 
understand aspects of clustering. The steady-state cluster size distributions of the model 
were obtained by numerically solving the rate equations describing the evolution of clusters. 
We have obtained scaling forms for the cluster-size distributions and the average cluster 
size. Our results agree well with experiments where clustering or chain formation have been 
observed. 
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Figure Captions 



Table I: Variation of the Neel (tjv), Brownian (tb) and effective (r e ) relaxation times as 
a function of the core radius r c of the magnetic nanoparticle. t# has also been evaluated 
for three values of the surfactant coating 5. The corresponding r e for these values is also 
evaluated. 



Table II: Calculated value of the ratio n and the corresponding estimated va 



cluster size for a variety of MNP used in experiments by Eberbeck, et al 
experimental parameters associated with each sample are also specified. 
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Figure 1: Variation of the cross-over radius r* as a function of the anisotropy constant K for 
several biologically and commercially used magnetic nanoparticles. The surfactant coating 
S and the temperature T were taken to be 2 nm and 300 K respectively in these evaluations. 

Figure 2: Variation of vs - 00 f° r monodisperse (open circles) and polydisperse (filled 

circles) manghemite samples when (a) r c < r* and Neel relaxation dominates, (b) f c « r* 
and both Neel and Brownian relaxation contribute and (c) f c > r* and Brownian relaxation 
dominates. The values of r c in (a), (b) and (c) have been chosen as 4, 8 and 12 nm 
respectively. The surfactant coating 5 and the temperature T were taken to be 2 nm and 
300 K respectively in these evaluations. 

Figure 3: Cole-Cole plots (open circles) for data of figure 2 when (a) r c < r* and Neel 
relaxation dominates, (b) r c ~ r* and both Neel and Brownian relaxation contribute 
and (c) r c > r* and Brownian relaxation dominates. The plots exhibit a two-humped 
structure when both relaxation times contribute (case b). The j3 and the r c values obtained 
from the Cole-Davidson equations (refer text in section II C) are also indicated in the figures. 

Figure 4: Variation of a as a function of /3 for three different anisotropy constants on a 
semi-logarithmic scale. 
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Figure 5: The reconstructed particle size distribution along with the particle size distri- 
bution obtained using TEM for comparisons (refer text in section II D for details). 



Figure 6: (a) Variation of the mean cluster size < k(R, t) > as a function of time t for ji = 
2 and different values of the ratio R specified in the figure. The inset shows the steady-state 
value of the mean cluster size < k(R) > as a function of R on a double logarithmic scale. 
The best fit line to data has a slope 0.85 ± 0.02 respectively, (b) Scaled data of (a) with 
the scaled indicated in the figure. The initial growth of clusters obeys a power 

law characterized by the exponent z = 0.38 ± 0.02. Similar data corresponding to /j, — 1 
is also indicated in the inset. For this data, the growth exponent z = 0.78±0.02 as indicated. 

Figure 7: (a) Steady-state distribution P(k). vs k for jj, — 2 for indicated values of ratio 
R. (b) Scaled data of (a) with the scaled clXIS cLS indicated in the figure. 



Figure 8: Data from reference [13J on average chain length as a function of time for different 
volume fractions replotted on scaled cLXGS cLS shown in the figure. The solid line corresponds 
to a slope of 0.78 obtained from a numerical solution of the aggregation- fragmentation model 
with /i = 1 (refer inset of figure 6b and corresponding text). 
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Table I 



r c [nm) 


tn (s) 

( s ) 
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Table II 



MNP Sample 


Core 


Shell 


Q 


< k > 


DDN128 


Fe 3 4 


carboxydextran 


0.5364 


2.3 


FluidMagD5 


Fe 3 4 


starch 


0.36 


2.0 


MagBSA 


Fe 3 4 


BSA 


0.676 


2.5 


Resovist 


Fe 3 4 


Carboydextran 


0.51 


2.4 
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Figure 6b 
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